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We report a study of the suppressed B meson decay B~ — > DK~ followed by D — > K + tt~, 
where D indicates a D° or D° state. The two decay paths interfere and provide information on the 
CP-violating angle (j>3- We use a data sample containing 657 x 10 6 BB pairs recorded at the T(4S) 
resonance with the Belle detector at the KEKB asymmetric-energy e + e~ storage ring. We do not 
find significant evidence for the mode B~ — > DK~ , D — > K + tv~ , and set an upper limit of rs < 0.19, 
where rs is the mag nitude of the ratio of amplitudes \A(B~ -> D°K~)/A(B~ -> D°K~)\. The 
decay B~ — > Dir~ , D — » K + -k~ is also analyzed as a reference, for which we observe a signal 
with 6.60" significance, and measure the charge asymmetry Ad-h to be — 0.02+q ^(stat) ±0.04(syst). 
In addition, the ratio B(B~ -» D°K~)/B{B- -> is measured to be [6.77 ± 0.23(stat) ± 

0.30(syst)] x 1CT 2 . 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 



Precise measurements of the parameters of the stan- 
dard model are fundamentally important and may re- 
veal newphysics. The Cabibbo-Kobayashi-Maskawa ma- 
trix [l|, Q consists of weak interaction parameters for 
the quark sector, one of which is the CP-violating angle 
cf> 3 = arg {—VadVub* /VcdVcb*)- Several proposed meth- 
ods for measuring </>3 exploit the interference between 
B~ -> D°K~ and B~ -> D°K~, where D° and D° decay 
to common final states 0,13- The effects of CP violation 
could be enhanced if the final state is chosen so that the 
interfering amplitudes have comparable magnitudes p|. 
The decay B~ DK~ , D -> K+ir~ (D = D° or D°) is 
a particularly useful mode, in which the color- favored B 
decay followed by the doubly Cabibbo-suppresscd D de- 
cay interferes with the color-suppressed B decay followed 
by the Cabibbo- favored D decay (Fig.[TJ. Previous stud- 
ies of this decay mode have not found a significant signal 
yield @, 0]. The decay B~ -> Dtt~, D -> K+-R- has a 
similar event topology and is Cabibbo-enhanced relative 
to the corresponding DK~ mode. Therefore this mode 
is an ideal control sample, while its CP asymmetry is 
expected to be negligible. 

In this analysis, we measure the ratios of the above sup- 
pressed decays relative to the favored decays B~ — > Dh~, 
D — > K~tt + , where h = K or w. The same selection cri- 
teria are used for the suppressed decays and the favored 
decays whenever possible in order to cancel systematic 
uncertainties. In this paper, charge conjugate reactions 
are implied except where otherwise mentioned; we de- 
note the suppressed decays B~ — ► Dh~ , D — > K + iv~ 
as B~ — > D sup h~ , and the favored decays B~ — ► Dh~, 
D — > K~ir + as B~ — > Df av h~ . Furthermore, a K~ or 
7T _ that originates directly from a B~ is referred to as 



B~ 



B~ 




FIG. 1: Diagrams for B~ -> DK~ , D -» K + n~ and B~ 
Dtt~ , D — > K + n~ decays. 

the "prompt" particle. 

The results are based on a data sample that contains 
657 x 10 6 BB pairs, collected with the Belle detector at 
the KEKB asymmetric-energy e + e~ (3.5 GeV on 8 GeV) 
collider Q operating at the T(45) resonance. 

The Belle detector is a large-solid-anglc magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of aero- 
gel threshold Cherenkov counters (ACC) , a barrel- like ar- 
rangement of time-of-flight scintillation counters (TOF), 
and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals (ECL) located inside a superconducting solenoid 
coil that provides a 1.5 T magnetic field. An iron flux re- 
turn located outside of the coil is instrumented to detect 
Kl mesons and to identify muons (KLM). The detector 
is described in detail elsewhere @. Two inner detector 
configurations were used. A 2.0 cm beam pipe and a 3- 
layer silicon vertex detector were used for the first sample 
of 152 x 1Q 6 BB pairs, while a 1.5 cm beam pipe, a 4-laycr 
silicon detector, and a small-cell inner drift chamber were 



used to record the remaining 505 x 10 6 BB pairs [13]. 

Neutral D meson candidates are reconstructed from 
pairs of oppositely charged tracks. For each track, we 
apply a particle identification requirement based on a 
K/tt likelihood ratio P(K/tt) = Ck/{Ck + C-n), where 
Ck and C n arc kaon and pion likelihoods, respectively. 
The likelihoods arc determined by the information from 
the ACC and TOF and specific ionization measurements 
from the CDC. We use the requirements P(K/ir) > 0.4 
and P(K/tt) < 0.7 for the kaon and pion candidates, 
respectively. The efficiency to identify a kaon (pion) is 
94%, while the probability that a pion (kaon) is misiden- 
tified as a kaon (pion) is about 10%. The systematic 
error in the K/tt selection efficiency is less than 1% for 
both kaons and pions. The invariant mass of the Kir pair 
must be within ±3cr of the nominal D mass [0]: 1.850 
GeV/c 2 < M(Kir) < 1.880 GeV/c 2 . To improve the 
momentum determinations, tracks from the D candidate 
are refitted with their invariant mass constrained to the 
nominal D mass. 

B meson candidates are reconstructed by combining 
a D candidate with a prompt charged hadron candi- 
date, for which the particle identification requirement 
P(K/tt) > 0.6 [P(K/tt) < 0.2] is used for B~ -> DK~ 
(B~ — > Dtt~). With this requirement, the efficiency to 
identify a kaon (pion) is 86% (81%), while the prob- 
ability that a pion (kaon) is misidentified as a kaon 
(pion) is about 5% (10%). The signal is identified by 
two kinematic variables, the energy difference AE = 
Eo+E h - — -Ebcam and the beam-energy-constrained mass 
M hc = V^-bcam _ \pd +Ph-\ 2 , where -E boam is the beam 
energy in the T(45) center-of-mass (cm.) frame. We re- 
quire Mbc to be within ±3c of the nominal B mass fllj ; 
namely, 5.271 GeV/c 2 < M bc < 5.287 GcV/c 2 . We 
then fit the AE distribution to extract the signal yield. 
In the rare cases where there is more than one candi- 
date in an event (0.3% for B~ — > D sup K~ and 0.7% for 
B~ — > D sup n~), we select the best candidate on the ba- 
sis of a x 2 determined from the difference between the 
measured and nominal values of M(Kir) and Mb c . 

The large background from the two jetlikc e + e~ — > qq 
(q = u,d,s,c) continuum processes is suppressed using 
variables that characterize the event topology. A Fisher 
discriminant [l2T ] made up of modified Fox- Wolfram mo- 
ments called the Super-Fox- Wolfram (SFW) [l3[ and 
cos 9b, where b is the angle of the B flight direction 
with respect to the beam axis in the cm. system, are 
employed. These two independent variables, SFW and 
cos9b, are combined to form likelihoods for signal (£ s ig) 
and for continuum background (£ C ont); we then construct 
a likelihood ratio TZ = £ s i g /(£ s i g + £ C ont)- We optimize 
the TZ requirement by maximizing S/y/ S + B, where S 
and B denote the expected numbers of signal and back- 
ground events in the signal region, using Monte Carlo 
samples. To estimate S, we consider only the contri- 
bution from B~ — ► D°K~ followed by D° — ► A'+7r~, 
where the value of tb of Eq. j6]) is taken to be 0.1. For 
B~ -> D aup K~ (B~ -> D sup n-) we require TZ > 0.90 



(TZ > 0.74), which retains 45% (70%) of the signal events 
and removes 99% (96%) of the continuum background. 
A similar TZ requirement is obtained if the optimization 
uses S/VB instead of S/y/S + B. 

For B~ — > D snp K~ , a possible background comes from 
B~ — > Dtt~ , D — > K + K~ , which has the same final state 
and the same position of the AE peak as the signal. We 
veto events that satisfy 1.840 GeV/c 2 < M{KK) < 1.890 
GeV/c 2 . After this veto, the estimated number of events 
that contribute to the signal yield is 0.22 ± 0.19. The 
favored decay B~ — > Df av /i _ can also produce a peaking 
background for the suppressed decay modes if both the 
kaon and the pion from the -Df av decay are misidentified 
and the particle assignments are interchanged. In or- 
der to remove this background, we veto events for which 
the invariant mass of the Ktt pair is inside the 1.865 
GeV/c 2 ± 0.020 GeV/c 2 window when the mass assign- 
ments are exchanged. After this requirement, we esti- 
mate that 0.17 ±0.13 (6.0 ±2.1) events contribute to the 
signal yield for B~ -> D sup K~ (B~ -> D sup 7r"). 

The signal yields are extracted using extended un- 
binned maximum likelihood fits to the AE distributions. 
For the signal, we use a sum of two Gaussians, where the 
parameters are determined by a fit to B~ — > _Df av 7r~. 
The same probability density function (PDF) is used for 
the signal peaks in all other modes; the validity of this 
assumption is verified by Monte Carlo studies. 

Backgrounds from B — > XK~ (X ^ Asup(fav)), such 
as B~ — > D*K~, can populate the negative AE region of 
the B~ — ► -D S up(fav)i^~ sample. The PDF for these back- 
grounds is obtained from the BB Monte Carlo samples, 
in which all known B and B meson decays are allowed. 
Similarly, backgrounds from B — > Xit~ (X ^ -D S up(fav))> 
such as B~ — > D*tt~ and B~ Dp~, can populate the 
negative AE region of the B~ — > -D SU p(fav) 7r_ sample, as 
well as the negative AE region of the B~ — > -D S up(fav)^ _ 
sample if the prompt pion is misidentified as a kaon. In 
the fit to B~ — > £'sup(fav) 7r_ the PDF of these back- 
grounds is obtained from the BB Monte Carlo sam- 
ples, while in the fit to B~ — > -D sup (fav)A'~ the PDF 
is obtained from data by assigning the kaon mass to the 
prompt pion track in the B~ — > -D S up(fav) 7r_ sample. The 
good quality of the fit of the B~ — ► Df av K~ data sample 
indicates the validity of this technique. 

The feed-across from the B~ — ► -D SU p(fav) 7r_ signal 
peak also appears in the fit to B~ — * D sup (f ay }K~ , where 
the prompt pion is misidentified as the kaon. The PDF 
is fixed from the fit to the B~ — > Df av 7r _ data sam- 
ple where the kaon mass is assigned to the prompt pion 
track. The shift caused by the incorrect mass assign- 
ment makes the shape of the AE distribution asymmet- 
ric, and thus we model the misidentification background 
as a sum of two asymmetric Gaussians, for which the left 
and the right sides have different widths. In the fit to 
B~ — > D sup K~, we fix the yields for the contributions 
from the B — > Xn~ background and the feed-across from 
the B~ ~ > D sup 7r~ signal peak, using the measured yields 



in the B~ — ► D sup iT~ sample scaled by the ratio of the 
B~ — > D[ av ir~ yields obtained in the — > D{ av K~ and 
£?~ — > Df av n~ analyses. 

The continuum background populates the entire A_E 
region, for which we use a linear function. The fit results 
are shown in Fig. [2l 



We obtain 



Rdk = [7.8±f 7 (stat)± 2 ;°(syst)] x 1(T 3 , 
R Dn = [3.40i°;i(stat)i°ii(syst)] X 10 



(2) 
(3) 



where the systematic errors (Table [TTJ) are subdivided as 
follows. 
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FIG. 2: AE distributions for (a) B~ -> 1 
B~ -> D sup ir-, (c) B~ -> D fm K~, and (d) B~ 
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Charge conjugate decays are included. In these plots, B~ — > 
DK~ components are shown by thicker dashed curves, and 
B~ — > Dir~ components are shown by thinner dashed curves. 
Backgrounds are shown by thicker dash-dotted curves (for 
B — > XK~), thinner dash-dotted curves (for B — » Xir~), 
and dotted curves (for the continuum). The sum of all com- 
ponents is shown by the solid curves. 



The 

K+TT- 

D sup K 



charmless decay B~ — > K + K~n~ (B~ — > 
7r~) can peak inside the signal region for B~ — > 
(£>~ — > Z? sup 7r~). For this background, we fit 
the AE distribution of events in the D mass sideband, 
defined as 0.020 GeV/c 2 < \M(Ktt) - 1.865 GeV/c 2 | < 
0.080 GeV/c 2 , and obtain an expected yield of -2.3 + 2.4 
(2.5 ± 4.5) events. We do not subtract this charmless 
contribution and instead include the uncertainties, +2.4 
(+4.5), in the systematic error. 

The signal yields (N Dh -) and the reconstruction ef- 
ficiencies (ejj/j-) for the decays B~ — > D sup h~ and 
B~ — > Df av h~ are listed in Table |TJ From the results, 
we calculate ratios of branching fractions, defined as 



Rdh 



B{B- 



) 



Nr 



B(B- 



D fa , v h ) N DtKrh -/e D{mh - 



(D 



(i) Fit: The uncertainties due to the PDFs of the B~ — > 

D sup (f av )/i _ decays and the qq background are ob- 
tained by varying the shape parameters by ±ler. 
Those due to the PDFs and yields of the back- 
grounds from B — ► XK~ and B — > Xtt~ are es- 
timated by fitting the AE distribution in the re- 
gion -0.05 GeV < AE < 0.15 GeV without in- 
cluding those contributions. The total fit error is 
the quadratic sum and 26% (3.1%) for Rdk {Rd-k)- 

(ii) Peaking backgrounds: The uncertainties due to the 

backgrounds which peak under the signal were 
described earlier, and the corresponding system- 
atic error in Rdk {Rd-k) is estimated to be i 2 5% 
(^5 3%)- This uncertainty is asymmetric because 
the uncertainty of the charmless background is 
taken only for the negative side. 

(iii) Efficiency: Monte Carlo statistics and the uncer- 
tainties in the efficiencies of particle identification 
requirements dominate the systematic error in de- 
tection efficiency, which is estimated to be 2.7% 
(2.5%) for Rdk (Rd v )- 

The total systematic error is the sum in quadrature of 
the above uncertainties. The possible fit bias is checked 
using a large number of pseudoexperiments and found to 
be negligible. 

The significances are estimated as In (Co/£ max ), 
where £ max is the maximum likelihood and Cq is the like- 
lihood when the signal yield is constrained to be zero. 
The distribution of the likelihood C is obtained by con- 
voluting the likelihood in the AE fit and an asymmetric 
Gaussian whose widths are the negative and positive sys- 
tematic errors. The results are shown in Table |U 

Since the signal for B~ — > D snp K~ is not signif- 
icant, we set an upper limit at the 90% confidence 



level (C.L.), Rdk < 1.8 x 10' 
calculated according to 



Jo 



+ This limit, 
£-(Rdk)<IRdk 



plimit 

= 0.9 



J £{Rdk)<IRdk ■ 

Using the values of Rdh obtained above and the B~ — > 
Df av h~ branching fractions from Ref. [Til ] . we determine 
the branching fractions for B~ — ► D sup h~ from 



B{BT 



B(B~ -> D i&y h-) x R Dh . (4) 



The results are summarized in Table |TJ For the B — > 
D S upK~ branching fraction, we set an upper limit at the 
90% C.L., B(B~ -> D sup K~) < 2.8 x 10~ 7 . Our branch- 



ing fraction for B~ — > _D sup 7r~ is consistent with the 
value expected from measured branching fractions for B 
and D decays (llj . 



TABLE I: Summary of the fit results. For the B~ — > D aup h~ signal yield, the contribution of peaking backgrounds has been 
subtracted. The first two errors on the measured branching fractions are statistical and systematic, respectively, and the third 
is due to the uncertainty in the B~ — > Df av ft~ branching fraction used for normalization. The last column shows the partial 
rate asymmetries Adh as explained in the text. 



Mode 



Efficiency (%) Signal yield Significance Branching fraction [90% C.L. upper limit] 



A L 



7+7.7 

1 -7.0 



™i±0.l)xfcr' 
[2.8 x 10" 7 1 



TUTS" 



B- D sup K~ 



B' 
B' 



D snp n 
D Cav K~ 

£>f av 7T~ 



f5.4±0.3 

23.f±0.4 
15.1±0.3 
22.8±0.4 



93.8 
f220 
27202 



-15.2 
-14.6 
+41 
40 
+ 177 
-176 



1.3a 
6.6a 



(1.2 



(6.2 9 ™ o o . 2 4i±0.24)xf0- 7 



-O.f 



-0.02^ 



±0.4 



±0.04 



TABLE II: Summary of the systematic uncertainties for Rdh 
and ADh- 



Source Rdk Rdtt Adk Ad-k 

Fit ±26% ±3.1% ±0.40 ±0.04 
Peaking backgrounds 1 2 5% i§'§% 

Efficiency ±2.7% ±2.5% 

Detector asymmetry • • • ■ ■ ■ ±0.01 ±0.01 



to the fits. Possible bias due to charge asymmetry of 
the detector is estimated using the B~ — > Df av ir~ con- 
trol sample for which the expected asymmetry is small. 
The peaking backgrounds are subtracted assuming no 
CP asymmetries. An assumption of 30% CP asymmetry 
in the peaking background would lead to a shift of 0.02 
in Ad-k- 



The ratio Rdk is related to <f>3 by 



Rdk 
where [l4{ 



r B + r D + ^r B r D cos cos 6 



I'D 



A(B- -► D°K~) 



A(B~ -> D°K~) 
A(D° -> K+TT-) 



A(D C 



0.0578 ±0.0008, 



(5) 

(6) 
(7) 



and 5b and Sd are the strong phase differences between 
the two B and D decay amplitudes, respectively. Using 
the above result, we obtain a conservative upper limit 
on tb as follows. For a given Rdk and in the relevant 
parameter ranges, tb is the largest when cos <r> 3 cos 5 = 

— 1 and td is maximal. Thus, we take cos 03 cos S = 

— 1 and a +2a shift in rn, and obtain tb < 0.19 which 
corresponds to the 90% upper limit on Rdk- 

We also measure the partial rate asymmetry Adh i n 
the B T — > D sup h T decays, 



Anh = 



_ B(B~ -> D mp h~) - B(B J 



n h~\ 

- Ly sup"' 



B[B- D sup h~) + B{B+ D snp h+) 



(8) 



by fitting the B~ and B + candidates with the asymmetry 
as one of the fitting parameters. The fit results are shown 
in Fig. [3] and included in Table HI We obtain 



Ai 



0.02+^°(stat) ± 0.04(syst) 



(9) 



and no significant constraint on Adk ■ The systematic 
errors (Table [TTJ) arc dominated by the uncertainties due 
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FIG. 3: AE distributions for (a) B~ -> D aup K', (b) B + -> 
D sup K+, (c) B~ -> D sup %-, and (d) B+ -> D Eup n+ . The 
curves show the B^ — > DsupA'^ component (thicker dashed 
curves), the B^ — > DsupTr^ component (thinner dashed 
curves), and the background components (thicker dash-dotted 
curves for B — * XK^, thinner dash-dotted curves for B — » 
Xn^, and dotted curves for the continuum), as well as the 
overall fit (solid curves). 



We also report the ratio 

B{B- -> D°K-) = N D ^ K -/e D!mK - 
B{B-^DH-) N Df ^-/e Dt ^- 1 ) 

to be [6.77 ± 0.23 (st at) ± 0.30(syst)] x 10~ 2 from the fit 
to B~ — > DfgjvK - and B~ — > Z)f av 7r~ J which is about 3cr 
lower than the current world average [111 ] . The systematic 
error is due to the uncertainties in the yield extractions 
(3.1%) and uncertainties in efficiency estimations (1.9%). 
The latter is dominated by the uncertainty in particle 
identification efficiency for prompt hadrons. 

In summary, using 657 x 10 6 BB pairs collected with 
the Belle detector, we report studies of the suppressed 
decay B~ — > D aup h~ {h = K,ir). No significant signal 
is observed for B~ — > D sup K~ and we set a 90% C.L. 
upper limit on the ratio of B decay amplitudes, rs < 
0.19. This result is consistent with the measurement of 



tb in the Dalitz plot analysis of the decay B~ — > DK~ , 
D — > KgTr + Tr~ [l5|, [l6[. For — * D sup 7r~, we observe a 
signal with 6.6er significance. We also report the charge 
asymmetry for B T — > D sup ir T and the ratio B(B~ — > 
D°K~)/B(B~ — > D°7r~). These results improve and 
supersede our previous results @, [T3| ■ 
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